Doping and temperature dependence of the pseudogap and Fermi arcs in cuprates 
from rf-CDW short-range fluctuations in the context of the t-J model 



Matias Bejas, Guillermo Buzon, Andres Greco, and Adriana Foussats 
Facultad de Ciencias Exactas, Ingenieria y Agrimensura and Instituto de Fisica Rosario, 
Universidad Nacional de Rosario and Consejo Nacional de Investigaciones CienUficas y Tecnicas, 
Avenida Pellegrini 250-2000 Rosario- Argentina. 
(Dated: January 12, 2013) 

At mean- field level the t-J model shows a phase diagram with close analogies to the phase diagram 
of hole doped cuprates. An order parameter associated with the flux or d charge- density wave (d- 
CDW) phase competes and coexists with superconductivity at low doping showing characteristics 
identified with the observed pseudogap in underdoped cuprates. In addition, in the d-CDW state the 
Fermi surface is reconstructed toward pockets with low spectral weight in the outer part, resembling 
the arcs observed in angle-resolved photoemission spectroscopy experiments. However, the d-CDW 
requires broken translational symmetry, a fact that is not completely accepted. Including self- 
energy corrections beyond the mean, field we found that the self-energy can be written as two 
distinct contributions. One of these (called E/z nx ) dominates at low energy and originates from 
the scattering between carriers and d-CDW fluctuations in proximity to the d-CDW instability. 
The second contribution (called T>r\) dominates at large energy and originates from the scattering 
between charge fluctuations under the constraint of non double occupancy. In this paper it is shown 
that E fiux is responsible for the origin of low-energy features in the spectral function as a pseudogap 
and Fermi arcs. The obtained doping and temperature dependence of the pseudogap and Fermi arcs 
is similar to that observed in experiments. At low energy, E#a gives an additional contribution to 
the closure of the pseudogap. 

PACS numbers: 74.72.-h, 71.10.Fd, 74.25.Jb, 79.60.-i 



I. INTRODUCTION 

The origin of the pseudogap (PG) phase in cuprates 
is one of the most important and unresolved issues in 
solid-state physics. 1 Several experimental techniques are 
used for studying this phase, and its main characteris- 
tics remain unclear. For instance, in the superconducting 
state, some experiments are consistent with the existence 
of only one gap while others are in agreement with two 
order (competing) parameters. Two main scenarios, pre- 
formed pairs above T c and two competing gaps, dispute 
the explanation of the PG (see Refs. and 111). 

Angle-resolved photoemission spectroscopy (ARPES) 
is presently a valuable tool for such research. 4 Surpris- 
ingly, in underdoped cuprates, ARPES shows, in the nor- 
mal state below a characteristic temperature T*, Fermi 
arcs£~— (FAs) (centered along the zone diagonal) in- 
stead of the full Fermi surface (FS) predicted by stan- 
dard solid-state physics. Despite the general consensus 
about the existence of FAs, their main characteristics 
are controversial, even at an experimental level. For in- 
stance, while some experiments suggest that the arcs are 
disconnected^"— others claim that the arcs are associated 
with pockets Ji"— 

The number of theoretical studies about the PG and 
FAs is too large for a complete listing. Among them, 
early results where the PG was discussed in the frame- 
work of the Born approximation and the spin-polaron 
description for the t-J model should be mentioned^ In 
addition, recent progress on dynamical cluster approxi- 
mation (DCA) show the presence of a pseudogap 23 ! 24 and 



Fermi arcs 2 ^"— in the two-dimensional Hubbard model. 

Recently, Norman et al. (Ref. 128l ) have summa- 
rized some of the relevant models proposed for discussing 
ARPES experiments. These models are semiphenomeno- 
logical or phenomenological, and the proposed Green 
function G{k 1 uo) has the following form: 



G~ 1 (k,uj) =u - e k + iT -E(fc,cj) 



(1) 



where T is a lifetime broadening, e& is the bare electronic 
dispersion, and E(fc,cj) is the self-energy, which can be 
written as 



E(/c,cj) 



CO 



(2) 



In Eq. (2) the phenomenological pseudogap is as- 
sumed to be d wave; = ^(cosk x — cosk y ). 
is model dependent; for instance: (a) = — c/c+q, 
where Q = (ir, n), in the d charge-density wave (d-CDW) 
model, 29 (b) is the nearest-neighbor term of the tight 
binding dispersion in the model proposed by Yang, Rice, 
and Zhang (YRZ)^S and (c) = e& in the d-wave pre- 
formed pairs model j 28 i 31 Although these models represent 
different physical situations, the experimental distinction 
between them is a big challenge. 

In Ref. [H it was concluded that d-CDW and YRZ 
models lead to predictions that are not compatible with 
experiments. For instance, these models lead to FAs 
whose length is temperature independent and deviates 
from the underlying FS in contrast to the experiments. 
Finally, it was also concluded that experiments are better 
described in the framework of the d-wave pairs model. 
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Since the early studies on high-Tc superconductivity, 
the t-J model has been shown to be a basic model for 
describing the physics of cuprates*^ This model, which 
can be considered a strong-coupling version of the Hub- 
bard model, 33,34 contains (potentially) the main ingre- 
dients for describing cuprates, i.e., antiferromagnetism 
at zero doping, a metallic phase at finite doping, strong 
tendency to d-wave superconductivity and several can- 
didates for the PG phase at low doping. Whether all 
these phases may be unambiguously associated to those 
known in cuprates is the big challenge for the t-J model. 
The number of analytical and numerical techniques in- 
troduced for studying this model is too large to discuss 
here. 

One analytical approach for treating the model is the 
large- N expansion where the two spin components are 
extended to N and an expansion in powers of the small 
parameter 1/N is performed. The advantage of this ap- 
proach is that the results are not perturbative in any 
model parameter and occur in strong coupling. For per- 
forming the large- N expansion, treatments based on the 
slave boson 3 - 5 - and Hubbard operators 36 were developed. 
However, evaluating fluctuations above mean field, as 
required for calculating dynamical self-energy effects, is 
not straightforward^ 5 - On the basis of the path-integral 
representation for Hubbard operators — the large- TV ap- 
proach to the t-J model was implemented, yielding pre- 
vious results 35,36 in leading order. At mean-field level 
(N = oo ) the well-known flux phase 3 -^"— (FP) instabil- 
ity at low doping was also reobtained. 37 In the FP a 
charge-density wave coexists with orbital currents in a 
staggered pattern and has the same momentum depen- 
dence of the superconducting state (d wave), allowing 
the identification of the FP with the pseudogap. In addi- 
tion, the FP scenario possesses the main properties to be 
identified with the phenomenological d-CDW proposal^ 
It is important to mention that the relevance of the FP 
for the physical case N = 2, for instance, in the form 
of a phase with strong d-wave short-range order, is un- 
der dispute. While some exact diagonalization results 42 
show the presence of the d-CBW phase, DCA (Ref. IZ3I ) 
and strong-coupling diagram technique 44 do not show the 
static long-range formation of the d-CDW state. In spite 
of this discussion it is important to note that the pre- 
dicted mean-field phase diagram 4 ^ has close similarities 
to the phase diagram of hole-doped cuprates where the 
FP competes and coexists with superconductivity. 

Importantly, the method developed in Ref. 37 allows 
us to go beyond the mean field and to compute self-energy 
renormalizations. Here, following Refs. EHI and l46l it will 
be shown that the doping and temperature dependence 
of the PG and FAs can be discussed after including self- 
energy effects in proximity to the FP instability, showing 
that d-CDW model is not inconsistent with the notion of 
arcs. 



This paper is organized as follows. In Sec. II we sum- 
marize the basic formalism. We show that the self-energy 
can be written in terms of two distinct contributions, 
Ttfiux and X^a- The mean- field phase diagram is dis- 
cussed together with the main characteristics of the self- 
energy. In Sec. Ill we describe the origin of the FAs 
and show that they are triggered by X/z n:c . Section III 
A discusses the topology of the FAs. Sections III B and 

III C discuss the temperature and doping dependence of 
the FAs, respectively. In Sec. Ill D we present the main 
characteristics of Y,fi ux at finite temperature. In Sec. 

IV we discuss the inclussion of Yir\. Section V presents 
discussion and conclusions. 



II. BASIC FRAMEWORK 

The large- N mean- field solution of the t-J model yields 
a quasiparticle (QP) dispersion^ 7 - 

e k = -2 ^ + rJ^j [cos(k x ) + cos(k y )\ 

-4 t' - cos(k x ) cos(ky) - /i, (3) 

where 5 is the doping away from half- filling, t, t 1 \ and 
J are hopping between nearest-neighbor, next-nearest- 
neighbor, and the nearest-neighbor Heisenberg coupling, 
respectively. The contribution r to the mean-field band 
and the chemical potential \i must be obtained self- 
consistently from 

T = TT C0S (k x )n F (e k ) (4) 
s k 

and 

(l-tf) = -^I>FfaO, (5) 

s k 

where n F is the Fermi factor and N s the number of sites. 

Equations (3) (5) define a homogeneous Fermi liquid 
(HFL) phase that, as discussed in Sec. I, is unstable 
against a flux phase or <i-CDW state at low doping. 

Beyond the mean field the computation of fluctuations 
in 0(1/ N) leads to the following expression for the self- 
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energy. 

lmE(k,u) = T" - ^2 h a (k,q,uj- e k - q )h b (k, q, u - e k - q ) 

q,a,b 

xIm[D ab (q,uj - e k - q )][n F (-e k - q ) + n B (u - ek-q)] 

(6) 

where n# is the Bose factor and the six-component vector 
h a (k,q,i/) is 



3 
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2/i 



Jr 



cos 



fe- f) cos (I) +cos(V I) cos (I) 



- Jr cos (a^ - y) ; - Jr cos (k y - ; Jr sin (k x - y) ; Jr sin - y 1 ) j . 



(7) 



The physical information contained in the vector 
h a {k,q,v) is as follows. The first component (called SR) 
is mainly dominated by the usual charge channel, the sec- 
ond component (called 5X) corresponds to the nondouble- 
occupancy constraint, and the last four components are 
driven by J. For the case J = the vector h a reduces to 



a two-component vector. 

In Eq. (6) D a b is a 6 x 6 matrix that contains contri- 
butions from the six different channels and their mixing. 

D~^(q,iuj n ) = [D^^iuJn)}- 1 - U ab (q,iuj n ) (8) 

where 



/ S 2 /2(V- 



D ab(^ iuJ n) 



J/2)[cos(q x ) 
5/2 







<Qy)] S/2 
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(9) 



Ck 



N s 



n F (e k ) , (10) 



with 



g(k,q,iuj n ) 



£k- 



tk 



(11) 



where iu n is the bosonic Matsubara frequency. Hereafter, 
t' jt — —0.35 and J ft = 0.3, which are suitable param- 
eters for cuprates, are used. The lattice constant a of 
the square lattice and t are considered to be a length 
unit and energy unit, respectively. In Eq.(j9j), V is the 
nearest-neighbor Coulomb repulsion. The role of V is to 
exclude phase separation. We choose V = 2 J. 

The instability of the mean-field solution occurs when 
det[D a b(q,iuj n = 0)] = (Ref. 1371 ) . For the present 
parameters, at T = 0, the instability takes place at S = 
S c ~ 0.23 for q ~ (vr, vr). It is important to note that D a b 
enters explicitly in the self-energy expression beyond the 
mean field [Eq.(6)]; thus, E probes the proximity to the 
instability at low uj and for momenta k — q near the FS. 

Since D a b contains contributions from six different 
channels and their mixing, it is important to isolate the 
most relevant channel dominating E near the instability. 
The eigenvector with zero eigenvalue of Dab takes the 
form ^ (0,0,0,0,-1,1) which is the eigenvector associ- 
ated to the FP instability. 37 Projecting E(fc,o;) on the 
FP eigenvector the following self-energy contribution is 



obtained^ 

ImZ f i ux (k,u) = -^^-y 2 (q,k)lmxfiux(q,u - e k - q ) 



N > , 

[n F (-e k - q ) + n B (uJ - e k - q )} 



(12) 



which shows the explicit contribution of the flux suscep- 
tibility 



Xfiux(q,oj) = [2 J r 2 - IL(q,u)] 1 
where U(q,uj) is an electronic polarizability 

>> 1 - 2 / 7,\ i n F(ek+ q ) - n F (e k )] 



(13) 



II(g, iuj n 



€k+q — £k 



(14) 



calculated with a form factor j(q,k) = 2r[sm(k x — 
q x /2) — sm.(k y — q y /2)\. Since the instability takes 
place at (vt,7t) the form factor j(q,k) transforms into 
~ [cos(fcaj) — cos(fc 2/ )], which indicates the d-wave char- 
acter of the FP. Thus, the mode associated with the FP 
instability plays an important role in E(fc,cj) at low dop- 
ing near S c . 

In Fig. 1, disregarding superconductivity, the solid 
line shows the temperature Tpp, which indicates the 
onset of FP instability, i.e., when the static (iuj n = 0) 
flux susceptibility [Eq.(13)] diverges. At T = a phase 
transition occurs at the quantum critical point (QCP) 



4 



0.03 



H 0.02 



0.01 




0.03 



J/t = 0.3 
t'/t = -0.35 



0.2 
doping 8 



0.4 



FIG. 1. (Color online) Phase diagram of the t-J model in 
the leading order oil/N expansion where superconductivity 
was discarded. The instability (solid line), marked by Tfp, 
separates the homogeneous Fermi liquid state from the flux 
or d-CDW state and terminates at the quantum critical point 
at 5 C ~ 0.23 at T — 0. The inset shows the imaginary part 
of the flux susceptibility vs. uj for 5 = 0.10 and for different 
temperatures [q = (7r,7r) is the momentum where the insta- 
bility occurs]. Approaching T F p from above the flux mode 
becomes better defined, accumulates weight, and softens. At 
T — Tfp the flux mode reaches uj = 0, freezing the d-CDW 
phase. This flux mode contributes to the self-energy leading 
to a pseudogap and Fermi arcs features as discussed in text. 
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FIG. 2. (Color online) The imaginary part of the two self- 
energy contributions, T,fi ux and E^a (E = T,fi ux + Era)- 
(a) Imaginary part of E/j ua . at T = for several dopings 
approaching the QCP S c ~ 0.23 for the antinodal Fermi wave 
vector kp N . (b) The same as (a) for E^a- Inset (i) shows the 
quasiparticle weight Z vs. doping for both contributions and 
for kp N and kp . While Z is weakly independent of doping and 
isotropic on the Fermi surface for E^a, for E/z nx Z is strongly 
anisotropic on the FS and strongly doping dependent. Note 
that Zfi ux for kp N tends to zero approaching the QCP. In 
addition, at large doping, the relevant contribution is E^a- 
Inset (ii) shows, for 6 — 0.24 and k — kp N , the — IhiEra in 
an extended uj scale, showing that the energy scale in E^a is 
large and of the order of t. 



placed at the critical doping S c . At Tpp a flux- mode 
[Imxfiux(Q = (7r, 7r),a;)] reaches uj = 0, freezing the FP. 
In the inset in Fig. 1, we have plotted lmxfiux(Q = 
(-7T, 7r), cj) for J = 0.10 for several temperatures approach- 
ing Tpp ft rsj 0.018, showing that when T — >• Tpp, a low 
energy d-wave flux mode becomes soft and accumulates 
large spectral weight. We have used the small broadening 
rj/t = 0.01 in the analytic continuation (iuj n uj + irj). 

Figure 2(a) shows — ImEfl ux (k, uj) at T = for several 
dopings at the antinodal Fermi wave vector kp N . At 
large doping 5 = 0.40, — ImEfl ux is weak and behaves 
as ~ uj 2 at low energy. Approaching S c (S = 0.26 and 
S = 0.24), — ImSflux increases, and the behavior at low 
energy is nearly linear in uo and develops structures at low 
energy uj/t ~ 0.1—0.2. Results for the nodal Fermi vector 
kp are not shown because they are nearly negligible due 
to the d-wave character of the flux instability. Inset (i) 
in Fig. 2(b) shows the QP weight Z at kp N (Zf£ x ) and 
at kp (Zfi ux ). While Z^ x is strongly doping dependent 
and tends to zero approaching 5 C , Z^ lux ~ 1 shows that 
Ttfiux is also strongly anisotropic on the FS. 

E fiux is written in terms of the flux or d-CDW suscep- 
tibility Xfiux(Q, which shows explicitly the role played 
by the soft flux mode with momentum (tt, tt) (see inset in 
Fig. 1). Therefore, near the antinode the QP on the FS 
is strongly distorted, leading to FA effects, as shown in 
Sec. Ill In addition, it is easy to check that the most im- 
portant J contribution to E(fc,u;) enters only via T>fi ux - 



In E(fc,u;), there is another contribution that is nearly 
independent of J. This contribution belongs to the usual 
charge SR and nondouble occupancy SX channels, (the 
first and second components of h a [Eq.(7)]) and can be 
written a o 45 i 48 

ImEp A (/c,cj) = - ^2 Im[D RR (q,uj - e k - q )] 

^ S q 

+ 2 Q lm[D X R(q,uj - e k - q )] 

+ lm[D\x(q,u- e k - q )]} 

x [n F (-e k - q ) + n B (u - e k - q )} , (15) 

where Q = ^(e k - q + uo + 2/x). 

Figure 2(b) shows -ImE RA (k,cj) at T = for kp N 
and for the same dopings as in Fig. 2(a). For kp, 
Im ER A (k, uj) (not shown) is nearly indistinguishable from 
results at kp N , showing that Ep A is rather isotropic on 
the FS. In addition, Ep A behaves as ~ uj 2 at low uj, and 
contrary to T,fi ux , there is no energy scale at low energy. 
In inset (i) we show Zp\ x and Zp^ x . These results show 
that the doping dependence of Ep A is weaker than T>fi ux . 
Note that Zr\ — ^ when 5—^0. It is important to note 
that there are no structures in ImER A (k,cj) at low en- 
ergy, and the main contribution appears at large ener- 
gies of the order of t [see inset (ii)]. Note also the strong 
asymmetry shown by Ep A that arises from nondouble- 
occupancy effects 4£ 

In summary, (a) Y,fi ux is strongly J and doping depen- 
dent, is strongly anisotropic on the FS, and contributes 
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at low energy, and (b) is nearly J and doping inde- 
pendent, is strongly isotropic on the FS, and contributes 
at large energy. Thus, E(fc,u;) can be written as the ad- 
dition of two well-decoupled channels. 

Im£(fe,o;) = ImE^(M) +lmY,fi ux (k,u) (16) 

Using the Kramers-Kronig relations, ReE(/c,co>) can 
be determined from ImE(fc,o;) and the spectral function 
A(k,(j), computed as usual: 

A(k u \ = _I ImgfoaQ 

1 ' ; 7T [u - e k - ReS(fc, w)] 2 + ImS(fc, w) 2 

(17) 

Before concluding this section it is important to note 
that at mean- field level the d-CDW picture leads below 
7Vp, where the translational symmetry is broken, to four 
hole pockets with low spectral weight in the outer side 
resembling the FAs 42 However, as discussed in Ref. [28|, 
this picture has conflicting points when compared with 
some ARPES data. We will show in Sec. Ill that the in- 
clusion of self-energy effects in proximity to the FP insta- 
bility provides a possible scenario for describing several 
ARPES features. Therefore, although at mean- field level 
the instability to the static d-CDW occurs below Tpp, 
the PG and FA formation do not require the long-range 
d-CDW state, but they do require the enhancement of 
fluctuations due to proximity effects. Thus, we are al- 
ways located in a homogeneous state with the presence 
of d-CDW fluctuations. 



III. Y> flux AND FERMI ARCS 
A. Topology of Fermi arcs 
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FIG. 3. (Color online) (a) Intensity of the spectral function 
at uj = vs kxjky for S = 0.10 and T/t = 0.025, above but 
close to TFp/t rsj 0.018. (b) The same as (a) but taken from 
the experimental results of Ref. 28 for comparison. As in 
the experiment, (a) shows a well-defined Fermi arc whose end 
does not turn away from the underlying FS (solid line). 
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FIG. 4. (Color online) (a) Intensity of the spectral function 
at the FS vs the Fermi surface angle (j) from the antinode 
(0 = 0o) to the node (0 = 45o) for 5 = 0.10 and T/t = 0.025. 
(b) The same as (a) but taken from the experimental results 
of Ref. [H for comparison. 



Since T,fi ux dominates at low energy we study here 
the spectral functions calculated with this contribution. 
In Sec. IV we show that the inclusion of Yir\ does 
not change the main conclusion obtained in this section. 
Figure 3(a) shows for 5 = 0.10 and T/t = 0.025 (above 
but close to Tpp) the spectral function intensity at zero 
energy vs k x ,k y . A well-defined FA is obtained. Similar 
to the experiment 2 ^ [Fig. 3(b)], the end of the arc does 
not turn away from the underlying FS, and there is no 
strong suppression of the intensity at the hot spots. 

In Fig. 4(a) the intensity of the spectral function on 
the FS is plotted as a function of the FS angle (j) [defined 
in Fig. 3(b)] from the antinode ((j) = 0o) to the node (</> = 
45o). As in the experiment 28 [Fig. 4(b)], the intensity 
monotonically decreases approaching the antinode but 
remains finite. In Fig. 5(a) energy distribution curves 
(EDC) on the underlying FS are plotted. In agreement 
with the experiment 28 [Fig. 5(b)], near the node, there 
are well-defined QP peaks; approaching the antinode, the 
spectral functions lose intensity at uj = 0, become broad, 
and develop a PG-like feature. The presence of a PG- 
like feature near the antinode means that the arc plotted 



in Fig. 3(a) is not simply related to a decrease in the 
intensity from the node to the antinode but that the FS 
near the antinode is gapped. 

Note that the present PG-like feature is not related to a 
true gap as in other models. It is developed dynamically 
in proximity to the d-CDW instability in the presence of 
strong and short-range d-CDW fluctuations. 

In summary, the effects described in Figs. 35 arise 
from self-energy effects due to the coupling between QPs 
and the soft flux mode (see inset in Fig. 1) in proximity 
to the FP instability (solid line in Fig. 1). Since the flux 
mode occurs mainly with momentum (-zr, 7r), the QP near 
the antinode is distorted, leading to FAs. Note also that 
since <i-CDW fluctuations are of short-range character, 
in the present picture, the translational symmetry is not 
broken. 



B. Temperature dependence of the Fermi arcs 

The temperature dependence of the length of the FAs 
is puzzling. In spite of different views and interpreta- 
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FIG. 5. (a) Energy distribution curves on the underlying FS 
from the antinode (top) to the node (bottom) for S = 0.10 
and T/t = 0.025. Near the node, well-defined QP peaks are 
observed. Moving from the node to the antinode, the intensity 
at uj — decreases, the peak becomes broad, and a pseudogap 
feature is developed, (b) The same as (a) but taken from the 
experimental results of Ref. 28 for comparison. 



tions most reports agree on the fact that the observed 
FAs depend on temperature. While there are reports 
that claim that the length of FAs collapse to one isolated 
point-^ (nodal metal) at T = 0, others suggest a less- 
strong temperature dependence ] 21 1 50 In models discussed 
in Ref. [28| the temperature dependence of the length of 
the arcs arises after assuming a temperature dependence 
for Afc or for the lifetime broadening T. In this subsec- 
tion it is shown that the temperature dependence of the 
FAs emerges, in the framework of the present approach, 
without adjustable parameters, showing that the temper- 
ature dependence of the length of the arcs is entangled 
to their origin. 

Figure 6 shows the plot, for S = 0.10, of FA for different 
temperatures. Clearly, the length of the arcs decreases 
when temperature decreases. We note that the temper- 
ature dependence of the arcs seems to be weaker than in 
some experiments 6 but closer to other o 21 i 5Q (this point 
is further discussed in Sec. IV). Beyond a quantitative 
comparison, the because no phenomenological parameter 
is assumed to be temperature dependent in the present 
model, the results can be considered satisfactory. Fig- 
ure 7 plots the spectral function intensity on the FS for 
several temperatures normalized to the intensity at kp . 
Consistent with the picture in Fig. 6, with decreasing 
temperature, the intensity is more concentrated around 
the node. 

Finally, Fig. 8 shows EDC at kp N for the same tem- 
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FIG. 6. (Color online) Fermi arc for 5 = 0.10 for several 
temperatures: (a) T/t = 0.021, (b) T/t = 0.025, and (c) 
T/t = 0.050. Similar to experiments, when the temperature 
increases, the length of the arc increases. 
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FIG. 7. (Color online) (a) Intensity of the spectral function 
(normalized to the intensity at the node ) at the FS vs. <j> 
for S = 0.10 and for several temperatures. When the temper- 
ature decreases toward Tfp, the intensity is more and more 
concentrated around the node. 
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FIG. 8. (Color online) EDC at kp N for 5 = 0.10 and for 
the same temperatures as in Fig. 7. As in experiments when 
temperature increases, the PG feature fades out. Although 
the leading edge of the pseudogap partially closes, a filling is 
also observed, as in the experiments. 
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FIG. 9. (Color online) Fermi arc for T/t = 0.025 and for 
several dopings: (a) S = 0.05, (b) S = 0.10, and (c) 5 = 0.15. 
As in experiments when doping increases toward overdoped, 
the length of the arc increases. 



peratures as in Fig. 7. Although the PG feature par- 
tially closes^ with increasing temperature, a filling is 
also observed. 5 7 This feature is in agreement with ex- 
periments and in contrast to results from mean-field cal- 
culations where only a closure is expected. 

It is worth mentioning that while the arcs discussed 
here are dynamically generated, they necessarily occur 
at finite temperature. The present approach should be 
distinguished from other ones 51 53 where a phenomeno- 
logical fitted susceptibility without explicit temperature 
dependence is proposed. 



C. Doping dependence of the Fermi arcs 

It is well known that with increasing doping, the PG 
feature closes 54,55 and, simultaneously, the length of the 
arcs increases X For describing this behavior, models dis- 
cussed in Ref. [HI need to assume a phenomenological 
doping dependence for the PG. In this subsection we will 
show that arcs whose length increases with increasing 
doping can be naturally discussed in the present context. 

In Fig. 9 the FA is shown for several dopings and for a 
fixed temperature T/t = 0.025. With increasing doping, 
the length of the arcs increases, in agreement with exper- 
iments. Figure 10 shows EDC at kp N for several dopings. 
When doping increases, the PG-like feature closes and, 
simultaneously, the intensity increases at uj = 0. 

In summary, in Sees. Ill B and III C it is shown that 
with increasing doping and temperature the PG-like fea- 
ture and the FA fade out like in the experiments. The 
origin for this behavior is easy to understand: By increas- 
ing doping and temperature we leave out the instability 
line Tpp. Then, the flux mode is less efficient, self-energy 
effects become weaker, and the long FS is smoothly re- 
covered. It is important to note that from our approach 
T* must be distinguished from a true phase transition. 
Here at T* > Tpp, where the PG features vanish, there is 
not a phase transition but a smooth crossover. 56 Finally, 
note that if t = OAeV, the energy scale for the pseudogap 
and temperature is of the order of the experiment. 
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FIG. 10. (Color online) EDC at kp N for T/t = 0.025 and for 
the same dopings as in Fig. 9. Like in experiments when dop- 
ing increases, the PG feature washes out in a way consistent 
with the increment of the length of the arc reported in Fig. 
9. 



D. Main characteristics of T,fi ux 

For a complete discussion about the origin of the arcs 
we have investigated the main characteristics of Y,fi ux . 
Figure 11 shows -ImS flux at k$ and k$ N for T/t = 0.025 
and S = 0.10 [Fig. 11(a)], 6 = 0.25 [Fig. 11(b)], and 
S = 0.40 [Fig. 11(c)]. At fcj)f, -ImS flux is smaller than 
for kp N , leading to a well-defined and nearly no renor- 
malized QP peak in the nodal direction [Fig. 5(a)]. How- 
ever, the behavior at kp N is very different, especially at 
low doping. Instead of a minimum at uj = 0, — ImSfl ux 
shows a maximum clearly observed for 5 = 0.10 [Fig. 
11(a)]. This behavior, which is in contrast to the ex- 
pected results from the usual many-body physics ^ 57 ! 58 is 
the main reason for the PG and FA formation. With in- 
creasing doping, the maximum at uj ~ washes out, and 
for large doping, — ImSfl ux develops the expected mini- 
mum at uj = 0. [See results for S = 0.40 in Fig. 11(c)]. 



IV. INCLUSION OF Y, RX 

It was shown (Fig. 2) that the energy scale in 
is much larger (~ t) than the energy scale in Ytfi ux . Al- 
though this fact implies (as shown in Sec. Ill) that Tifi ux 
is the relevant contribution for triggering the low-energy 
PG features, in this section, for completeness, we discuss 
the role of in the spectral functions. It was discussed 
in Sec. IIIB that the PG closes and fills smoothly with 
increasing temperature (see Fig. 8). In this section we 
show that the only role of including Yir\ is to improve 
the vanishing of the PG. 59 

Figure 12(a) shows EDC for 5 = 0.10 for several tem- 
peratures at kp N . This figure shows that with increasing 
temperature the PG fills and a peak at uj = emerges 
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FIG. 11. (Color online) -ImS flux for T/t = 0.025 at k$ N and 
k% for (a) 8 = 0.10, (b) 5 = 0.25, and (c) 5 = 0.40. For all 
dopings, — ImSflux at kp (dashed line) is smaller than the cor- 
responding results at hp and shows the expected minimum 
at uj = 0. However, the behavior at kp N is very different. 
With decreasing doping, — ImEfl ux (solid line) increases, and 
a maximum, instead of a minimum, is developed at uj ~ 0. 
This behavior, which occurs only at finite temperature, is the 
cause for the dynamical generation of the arcs and the PG 
feature. For large doping (see results for S = 0.40), — ImEfl ux 
is small and depicts the expected behavior from the usual 
many-body theory, i.e., it has a minimum at uj = 0. This 
behavior is consistent with the fact that no arcs and no PG 
features are obtained for large doping. 



at T/t ~ 0.035. Note that in Fig. 8, where only E fiux 
was considered, even at the high temperature T/t = 0.1 
the maximum of A(k, uj) is not yet fully formed at uj = 0. 
In Fig. 12(c) we have reproduced, for comparison, the 
experimental results, 6 showing qualitative agreement be- 
tween theory and experiment. In Fig. 12(b) we plot, for 
T/t = 0.035, A(k,cj) for kp N (solid line) and k^ (dashed 
line). Althought the entire FS is ungapped at this tem- 
perature, the QP are better defined near the node, as in 
the experiment. 54 

Figure 13 shows the spectral function intensity at uj = 
vs k Xl k y for 5 — 0.10 and for the same temperatures 
as in Fig. 6. At low temperatures a FA is obtained, and 
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FIG. 12. (Color online) (a) EDC at k£ N for S = 0.10 and 
for several temperatures calculated using both contributions, 
T>r\ and S/^^. Note that the inclusion of T,r\ does not 
change the main conclusion obtained when only ^fi ux is con- 
sidered. As in the experiments, when temperature increases, 
the PG feature washes out. As in Fig. 8, although the leading 
edge of the pseudogap closes, a filling is also observed. Note 
that different from the calculation with only Yifi ux (Fig. 8), at 
T/t — 0.035 a full peak is recovered, (b) Spectral functions at 
kp N and k^ for T/t = 0.035. (c) Experimental results taken 
from Ref. [y for comparison. 
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FIG. 13. (Color online) Fermi arc for the same temperatures 
as in Fig. 6, (a) T/t = 0.021, (b) T/t = 0.025, and (c) 
T/t = 0.050, but calculated using both Hr\ and ^fi ux - As 
in Fig. 6, when temperature increases the length of the arc 
increases. 



its length increases with increasing T. Note that while 
the FS is expected to be gapped near the antinode for 
T/t = 0.025 (dot-dashed line in Fig. 12), for T/t = 0.050 
the full FS is ungapped. 

In summary, does not modify the main conclusion 
obtained in Sec. III. We have shown that its inclusion en- 
hances the pseudogap closing and filling, and contributes 
to a faster reconstruction of the entire FS with increasing 
temperature. 



V. DISCUSSION AND CONCLUSION 

The large- approach in the t-J model leads, beyond 
the mean-field level, to two distinct dynamical self-energy 
contributions, namely, and Y,fi ux . In this paper we 
have analyzed the role of these contributions in ARPES. 
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The main characteristics of Tifi ux are the following. 
^fiux is strongly anisotropic on the FS, strongly doping 
dependent, and dominated by J (if J = T,fi ux is neg- 
ligible), and it contributes at low energy. Thus, T,fi ux 
is the relevant contribution for describing the Fermi arcs 
and pseudogap features. 

The fact that E fi ux is mainly dominated by J may be 
understood as follows. At mean- field level the t-J model 
shows (and only for finite J) the flux or d-CDW phase 
below a temperature Tpp. Tpp decreases with increas- 
ing doping, approaching the QCP at S c and T = 0. In 
the proximity of Tpp, <i-CDW fluctuations enter T>fi ux 
[Eq.(12)]. Since d-CDW fluctuations favor scattering be- 
tween electrons with momentum transfer q ~ (tt, 7r), the 
FS near the antinode is gaped, leading to Fermi arcs be- 
ing dynamically generated. With increasing doping and 
temperature beyond S c and Tpp, respectively, d-CDW 
fluctuations become weak, and the Fermi arcs and the 
pseudogap wash out, in agreement with experiments. 

It is important to note that the present picture does 
not require any phenomenological parametrization for 
the pseudogap or the lifetime broadening and their tem- 
perature and doping dependence. It is only necessary 
to be in the proximity of the d-CDW instability or in a 
situation with strong short-range fluctuations. In other 
words, under the present approach Fermi arcs originate 
dynamically due to the interaction between carriers and 
short-range and short-living d-CDW fluctuations, imply- 
ing that long-range order is not broken. 

The present picture has similarities with some phe- 
nomenological approaches 53 where the pseudogap and 
Fermi arcs are described in a scenario where fermions 
interact with bosonic fluctuations of some special order. 
Importantly, our description offers a microscopic deriva- 
tion from the t-J model, and, as a corollary, the fluctu- 
ating spectrum is obtained with no assumptions of any 
fitted phenomenological parameter, such as coupling, cor- 
relation length, or bosonic frequency. Note that near the 
flux instability the flux mode (inset in Fig. 1) is over- 
damped and intrinsically temperature dependent and can 
not be easily considered as an Einstein mode as in other 
approaches^ - — 

A recent ARPES experiment^ suggests a similar sce- 
nario to that presented here, i.e., density wave fluctua- 
tions without long-range order. As in that experiment, 
in our theory, the existence (and persistence with de- 
creasing temperature) of broad and gapped spectral fea- 
tures near the antinode means that we are not sitting 
in a phase with long-range order. It is worth mention- 
ing that under the present approach, below the mean- 
field temperature Tpp the long-range d-CDW order oc- 
curs; that is, a true gap is formed, and sharp spectral 
features are expected with the corresponding reconstruc- 
tion of the FS in the form of pockets. 46 From an exper- 
imental point of view the existence of long-range order 
in underdoped cuprates is controversial and is tied to 
the following facts, (a) Some ARPES experiments show 
well-defined spectral peaks near the antinode in the su- 



perconducting state, while others show broad structures 
(see Ref. 60 and references therein), (b) While some 
experiments support the existence of a second order pa- 
rameter, distinct from but coexisting (and competing) 
with superconductivity^"— others claim to observe only 
one gap feature* 5 — 8 (c) While quantum oscillations 61 and 
some ARPES experiments show a reconstruction of the 
FS in the form of pockets ji 7 - - — other reports show only 
arcsA 7 -^ Although it is not our aim to solve these puz- 
zles (which requires more theoretical and experimental 
work) , we have shown that several aspects related to the 
Fermi arc phenomenology can be explained by d-CDW 
proximity effects, showing that this picture is not neces- 
sarily inconsistent with the notion of arcs. 

The characteristics of E#a are very different from those 
of Ytfiux. E^a is dominated by the usual charge chan- 
nel and (nearly) independent of J. Thus, this is the 
relevant contribution for the J = case. In addition, 
it is strongly asymmetric in uj around the FS due to 
nondouble-occupancy effects, rather isotropic on the FS, 
and rather constant as a function of doping (for low to 
intermediate doping). 48 Finally, it contributes at large 
energy of the order of t. Although E#a is not responsible 
for the pseudogap and Fermi arc formation, it gives an 
additional contribution to the temperature vanishing of 
the pseudogap. 

E^a and E fi ux may also play a role in describing other 
experiments in cuprates. (a) Since E^a shows high- en- 
ergy contributions, it leads, in the spectral functions, to 
incoherent structures at high binding energy, which of- 
fers a possible explanatio n 48 ! 62 for the high-energy fea- 
tures or waterfall effects observed in cuprates^ 3 -"— Other 
theoretica l 67 ! 68 and experimental 66 reports show a simi- 
lar conclusion, (b) The existence of two self-energy con- 
tributions is also consistent with recent angle-dependent 
magnetoresistance (ADMR) experiments^ - — These ex- 
periments show two different inelastic scattering rates 
with similar characteristics to the self-energy behavior 
discussed here, i.e., a strongly-doping-dependent and 
anisotropic scattering rate on the FS and another one 
that is weakly doping dependent and isotropic on the 
FS. Recently, ADMR experiments were discussed in the 
context of the present approach* 7 ^ 

Here we want to comment on the recent progress on 
DCA. As discussed in Sec. I DCA shows the presence of 
a pseudogap 2 — *2i and Fermi arcs^"— We wish to men- 
tion here the similarities between our results and those in 
DCA. For instance, the pole feature at uj ~ and near the 
antinode that occurs in JmE (Fig. 11), which diminishes 
with increasing temperature and doping, is in remarkable 
agreement with similar results discussed in Ref. |27]. This 
behavior for the self-energy leads also to a similar doping 
and temperature dependence for the PG and FAs. Note 
that in Ref. [53 the temperature filling of the PG as dis- 
cussed in the present paperwas also obtained. We note 
again that our results do not require the static long-range 
order of the d-CDW. What is needed is the enhancement 
of the d-CDW susceptibility due to fluctuations, as can 
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be seen in the inset of Fig. 1. Interestingly, although 
the static d-CDW state was not found in Ref. |43|, an 
enhancement of the d-CDW susceptibility was obtained. 
Finally, it is worth mentioning that although the origin 
of the PG and FAs is presumably of antiferromagnetic 
nature ^ 23 ! 24 a recent report is not conclusive about this 
affirmation. It is the aim of the present paper to show 
that d-CDW fluctuations may contribute to the PG and 
FAs formation. 

Although of one could certainly suppose that the large- 
N is a particular approximation and some results may 
depend on its details, we think that our theory contains 
features that can be expected, qualitatively, in cuprates 
and in the t-J model. Since the low-energy pseudogap 
feature increases with decreasing doping, it is reasonable 
to think that the pseudogap is associated with the same 
energy scale as the antiferromagnetism, i.e., J. This fact 



is contained in Y*fi ux . On the other hand, there is a 
larger energy scale, the hopping t, which, together with 
nondouble-occupancy effects, enters through Yir\. 

Finally, it is worth mentioning that besides d-CDW, 
other instabilities like stripes, 73 antiferromagnetism^i 
and Pomeranchuk 75 have been proposed to exist at low 
doping in cuprates. Thus, it is important to perform 
similar calculations for those cases and compare different 
predictions. 
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